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The centrality dependence of thermal parameters, characterizing the hadron multiplicities, is de- 
termined phenomenologically for lead-on-lead collisions at CERN-SPS for a beam energy of 158 
AGeV. The strangeness equilibration factor shows a clear, approximately linear, increase with in- 
creasing centrality, while the freeze-out temperature and chemical potential remain constant. 
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It has been shown that the abundances of different hadronic species in the final state of relativistic heavy-ion 
collisions can be well described by statistical-thermal models (cf. 0, ||, |^ and further references therein). In such a 
way, the observed multiplicities of a large number of hadrons can be reproduced by a small number of parameters 
J — 1^7^ I such as the temperature, the baryon chemical potential, the volume, and a parameter Q measuring the degree 
T-H , of equilibration of strange particles. Such a description can be justified for multiplicities measured over the whole 
phase-space, since many dynamical effects cancel out in ratios of hadron yields; in particular, effects due to flow can 
' be shown to disappear if the freeze-out surface is characterized by a single temperature and chemical potential |j . 
^ , As a matter of fact we mention that electromagnetic signals, i.e. real and virtual photons, observed in central 
' heavy-ion collisions at CERN Super Proton Synchrotron (SPS) energies can also be described astonishingly well 
, within the framework of a thermal fireball model . This is not necessarily a proof for (global) thermalization [|j , 
' but shows that the use of thermal models allows a very economical description of a large set of observables by a few 
T-H ' characteristic parameters. 

^-H . It is the subject of the present paper to pursue this idea and to analyze the centrality dependence of the thermal 
' parameters describing hadron multiplicities. This will provide further information about the effects of the size of the 
, excited strongly interacting system and help in the systematic understanding of the experimental data. We will show 
' that the thermal model is able to describe the available data for various centrality classes at one beam energy. 
^ In order to have a sound basis for the application of the thermal model we rely as much as possible on fully 
^ • integrated particle multiplicities. For this reason we concentrate our efforts on the analysis of results obtained by the 
^ , NA49 collaboration using centrality selected fixed-target Pb + Pb collisions at a beam energy of 158 GeV per 
nucleon, which are here analyzed within the framework of the thermal model § . The data jl^ consists of centrality 
binned multiplicities of the following hadronic species: p, (tt) = (tt"*" -|-7r~)/2, if + and , as a function of the number 
of participants. Following |l^, |ll[, these centrality classes correspond to mean impact parameters (in units of fm) 2.2 
y (bin I), 4.3 (bin II), 6.0 (bin III), 7.4 (bin IV), 8.9 (bin V), and 10.7 (bin VI). The mean number of participants is 
listed in Table 1 . Notice that the number of hadronic species at our disposal is quite restricted in comparison with 
other available data sets for central collisions. We investigate here the general trend for changing centrality, relying 
on the data set given in [p^ . It should be emphasized that we have not included in our analysis the multiplicity 
of protons. At larger impact parameters many protons are spectators which should not be included in the thermal 
model analysis. 

For completeness we briefly recall the basic features of the statistical-thermal model used to extract the chemical 
freeze-out parameters. A more detailed description can be found in p^ . 

• An ideal gas of all hadrons and hadronic resonances listed by the particle data group |T3| was used. 



The grand-canonical ensemble description was used throughout with the additional parameter 7s (cf. Q) in- 
cluded to account for the incomplete equilibration in the strange sector. 

The correct quantum statistics was used for each particle. 

Resonances were described with a Breit-Wigner distribution over an interval [m - Sm, m + 2T] where 
dm — min[mthicshoid, 2r]. The resonance widths F were taken from pM. 



• All resonances were allowed to decay with branching ratios taken from |13t . Weak decays were excluded since 
the data in were corrected accordingly. 
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• No excluded- volume corrections (cf. [Q) were incorporated. This does not affect the intensive thermal param- 
eters (T, /is, etc.), however, corrections to the densities will be called for, in particular the baryon and the 
energy densities. 

According to our proposition, at chemical freeze-out the system is completely specified by its temperature T, 
effective volume Vcs, strangeness suppression factor 7^ and chemical potentials ^b, Ms ^^id /ig. The initial conditions 
of the collision fix and /ig for a given ^b', Ms is fixed by the requirement of strangeness conservation, while /ig is 
determined by the constraint that the ratio of the total charge to total baryon number of the system is conserved. 

The chi-square minimization routine MINUIT was used to fit the chemical freeze-out parameters using 
= X]j Br{j i) Tj^'P''""-' v^jjT where Br{i —> i) is the branching ratio of the decay of hadron species j into species i, 
and n^''"^'™'' denotes the primary density. The values of T, /is, 7^ and VcS were used as free fit parameters adjusted 
to the fully integrated yields p, , , (tt) and Npart from jl^. The results obtained for each of the six centrality 
bins are listed in Table 1 and displayed in Figs. 1-4. The number of participants, TVpart, is identified with the baryon 
number of the hadron gas. 

As seen in Table 1, good agreement exists between the trends in the data and our model results. As the number of 
pions has a very small experimental error, the minimization routine zooms in on them and fixes the parameters so as 
to reproduce this number as accurately as possible. On the contrary, the number of participants, Npart, has a large 
error and the minimization routine gives much less weight to them, hence the agreement is not as good. 

Confident that the thermal model reproduces the available hadron multiplicities in each of the centrality bins with 
sufficient accuracy, let us consider the trends of the thermal model parameters. The chemical freeze-out temperature 
and the chemical potential stay fairly constant for the various centralities, as seen in Figs. 1 and 2. 

A comparison of T and /i b with the values obtained in the recent analysis of |^ for the most central collisions shows 
that they are consistent with the values of the present analysis. The results from |^] are shown as a band in Figs. 1 
and 2. As proven in ||^, the inclusion or omission of certain hadron species can change considerably the extracted 
values of T and /i^. We stress however that our analysis, due to the restricted available data, focuses on the trends 
with changing centrality. The temperature, T, and the baryon chemical potential, /is, do not show any noticeable 
dependence on centrality. 

There arc two thermal parameters which exhibit a pronounced dependence on centrality: the strangeness equili- 
bration factor 7s increases approximately linearly and then saturates with increasing centrality, as shown in Fig. 3; 
similarly, the radius R of the hadron gas, which can be readily compared to the radius of a static Pb nucleus increases 
approximately linearly with the number of participants, Npart, as seen in Fig 4. 

Fig. 5 shows the centrality dependence of various thermodynamic state variables, such as the energy per hadron 
{E)/{N), the energy density, the baryon density, and the entropy per baryon S/B. They all remain fairly independent 
of the centrality. In contrast, the Wroblewski factor defined as 

A, = -4^4^ (0.1) 

{uu) + {dd) ^ ' 

increases nearly linearly with increasing centrality. 

The predictions of other hadron multiplicities, such as K^, (p, S^, and A, A can serve as a further test of the 
thermal model. These are not yet at our disposal in such a centrality binned way as the yields of tt*, K^, and p but, 
may become accessible in further data analyses. We therefore present our prediction in Table ||. 

In summary, the analysis of the thermal parameters, describing the integrated yields of tt*, , and p as obtained 
by the NA49 experiment [|lO|, shows that the radius of the fireball increases linearly with increasing centrality. Also, 
the strangeness parameter 7^ increases, i.e., strange particle multiplicities approach chemical equilibrium. In contrast, 
the temperature and the baryon chemical potential do not change with centrality. No in-medium modifications are 
needed to describe the above quoted hadron yields. 
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TABLE L Comparison of experimental hadron yields with the results of the thermal model with parameters as displayed in 
Figs. 1 - 4. 
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TABLE II: Prediction of further hadron multiplicities using the parameters as displayed in Figs. 1-4. 
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FIG. 1: Centrality dependence of the chemical freeze-out temperature T as a function of the mean participant number. The 
result of [H] is indicated by the band: T = 158.1 ± 3.2 MeV. 
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FIG. 2: Centrality dependence of the baryon chemical potential fj,B- The result of is indicated by the band: fiB = 238 ± 13 
MeV. 
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FIG. 3: Centrality dependence of the strangeness equilibration factor 7s 
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FIG. 4: Centrality dependence of the fireball radius R. 
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FIG. 5: Centrality dependence of various thermodynamic state variables. The entropy per baryon is scaled by a factor 1/20. 
The baryon density is in units of 1/fm^, the energy density is in units of GeV/fm^ and are both uncorrected for excluded 
volume effects. The energy per hadron }E) / (N) is in units of GeV. 



